Epigenetics & Chromatin

Nuclear organization is dynamic and originates from the fundamental subunit of chromatin, the nucleosome. Post-translational modification of nucleosomal histones, particularly within intrinsically disordered histone tail regions, provides a dynamic regulatory mechanism of accessibility for chromatin-templated processes. While the epigenomic impacts of lysine acetylation and serine phosphorylation in the histone H3 tail are well-known, how these charge-altering post-translational modifications (PTMs) alter nucleosomal tail conformational dynamics remains incompletely characterized. Given that the functional implications of these PTMs are, at least in part, a consequence of modified nucleosome conformation, systematically cataloging the impact of histone PTMs on nucleosome dynamics provides crucial insight into both baseline cellular activity and epigenetic dysregulation that occurs in disease. Previously, our lab demonstrated that arginine citrullination mimetics lead to regional increases in H3 tail dynamics within nucleosome core particles. Here, we performed nuclear magnetic resonance spin relaxation experiments to investigate the effects of lysine acetylation and serine phosphorylation on H3 tail picosecond-nanosecond (ps-ns) dynamics. Using lysine-to-glutamine and serine-to-glutamate mutations as acetyllysine and phosphoserine mimetics, respectively, we found that these PTMs increase ps-ns conformational dynamics regionally around the PTM site, with a position-dependent effect. Additionally, we show that the type of PTM influences the extent of these increases: in general, the effect of mimetics trends in the order of phosphorylation [&le;] acetylation < citrullination, suggesting a tunable method for altering histone tail dynamics. Taken together, these results illustrate the role of nucleosome conformational dynamics in conveying the effects of epigenomic PTMs, elucidating a mechanism of the histone language.

Long-read nanopore sequencing enables simultaneous detection of germline variation and native DNA base modifications on individual DNA molecules, providing a unique opportunity to investigate allele-specific epigenetic regulation. Here, we performed whole-genome nanopore sequencing on normal and tumor prostate tissues to characterize differential methylation, methylation entropy, and allele-specific methylation (ASM) associated with noncoding genetic variants. Genome-wide analysis identified extensive cancer-associated differentially methylated regions (DMRs), with hypermethylated DMRs significantly enriched near transcription start sites and transcriptional regulatory regions. Integration with transcriptomic datasets revealed strong inverse relationships between promoter methylation and gene expression, while 5-hydroxymethylcytosine (5hmC) levels positively correlated with transcriptional activity across gene bodies. Using fragment-level methylation patterns enabled by long-read sequencing, we further quantified methylation entropy incorporating both 5mCG and 5hmCG states. Cancer-hypermethylated DMRs exhibited markedly reduced entropy, consistent with clonal fixation of methylation states during tumor progression. Entropy profiling across chromatin annotations demonstrated maximal epigenetic heterogeneity at partially modified enhancer-associated regions. To investigate cis-regulatory genetic effects, we developed a simple ASM framework (nanoASM) that can partition sequencing reads by allelic state and identifies allele-specific DMRs directly from long-read data. Compared with conventional population-level mQTL analysis, ASM demonstrated substantially improved statistical efficiency by leveraging within-individual contrasts and reducing sample-level heterogeneity. Although germline single nucleotide polymorphisms (SNPs) were largely shared between normal and tumor tissues, ASM patterns differed substantially, with tumor-associated ASM regions displaying significantly larger genomic span and stronger allelic methylation differences. Comparative analysis with TCGA prostate mQTL and GTEx prostate eQTL datasets demonstrated substantial concordance between ASM directionality and downstream transcriptional effects, particularly for variants located within DMRs and near transcription start sites. At the IRX4 prostate cancer risk locus, ASM identified an androgen-responsive regulatory domain overlapping AR ChIP-seq and H3K27ac peaks, nominating rs6885084 as a candidate functional variant. At the PSCA locus, ASM anchored by rs4736369 was associated with allele-specific methylation, chromatin activation, transcript abundance, and isoform usage. Together, these findings establish nanopore-based ASM analysis as a powerful approach for resolving functional noncoding variants and their regulatory domains they control in prostate cancer.

Nucleosomes are fundamental units of DNA packaging and gene regulation in eukaryotes. In mammalian sperm, most nucleosomes are replaced by protamines causing extreme chromatin compaction. Various epigenomic studies reported conflicting results on the distribution of residual nucleosomes in mammalian sperm, questioning their potential role in mediating intergenerational inheritance of paternal epigenetic information. Here we performed single-molecule footprinting through Nucleosome Occupancy and Methylome (NOMe) sequencing and applied the Bayesian statistical model nomeR to determine frequencies of nucleosome removal and retention at 103 specific genomic regions in thousands of developing haploid spermatids and mature spermatozoa of mice. While we readily detected footprints of nucleosomes and the transcription factor CTCF in round spermatids, chromatin became transiently highly accessible in elongating spermatids with loss of such footprints, indicating extensive chromatin reprogramming during spermiogenesis. In mature sperm, following nuclear decondensation with recombinant nucleoplasmin, we measured nucleosome occupancy frequencies ranging ~1.2 to 1.7% at mouse loci. In human sperm, nucleosome occupancy varied between ~2.3 to 4.5% at 163 genomic loci profiled. Contrasting mice, chromatin in ~25% of human sperm was accessible upon reducing disulfide bonds between protamines arguing for species specific protamine packaging. Our findings support a stochastic rather than programmed potential role of residual nucleosomes in mammalian sperm in regulating paternal gene expression during ensuing embryonic development.

The information content of the genome has been extensively analyzed. However, a comparable quantitative framework for DNA methylation is still lacking. Without such quantification, the magnitude of this regulatory and dynamic epigenetic structure remains conceptually imprecise, even though methylation dysregulation is strongly linked to disease-related phenotypes and altered cellular identity. Here we address this gap by applying Shannon information theory to DNA methylation. We first consider methylation marks as binary or probabilistic regulatory states and estimate the theoretical upper-bound information capacity of the human methylome under simplifying assumptions. We then progressively refine this estimate by incorporating biologically relevant constraints, including methylation bias, bimodal methylation distributions, local CpG correlation, genomic regulatory class, and cell-type-discriminative methylation patterns. This approach allows us to distinguish between theoretical methylation capacity, statistical methylation entropy, and biologically interpretable regulatory information. Finally, we consider methylation information from a discriminative perspective, analyzing its contribution to distinguishing cell types and regulatory cellular states. Within this framework, mutual information between methylation patterns and cell identity provides a biologically constrained estimate of methylations role as an epigenetic identity code. Our layered analysis reconciles megabit-scale methylome capacity with compact, biologically interpretable identity signatures. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=113 SRC="FIGDIR/small/735086v1_ufig1.gif" ALT="Figure 1"> View larger version (73K): org.highwire.dtl.DTLVardef@f0f0fdorg.highwire.dtl.DTLVardef@5d8a1eorg.highwire.dtl.DTLVardef@116debdorg.highwire.dtl.DTLVardef@79530e_HPS_FORMAT_FIGEXP M_FIG C_FIG

Chromatin-based genomics data are essential for understanding genome regulation and the mechanisms underlying epigenetic memory. Recent methods such as ChOR-seq and SCAR-seq assess histone modifications and chromatin-associated proteins during and after replication, capturing chromatin states that contribute to memory across cell divisions. Current tools for chromatin data analysis lack scalability and reproducibility across computing infrastructures, offer limited parameters, and are applicable only to a few sequencing techniques, ignoring the information from nascent chromatin assays. To address these challenges, we developed CREPAS, a Nextflow pipeline for analyzing nascent and parental chromatin sequencing data, including ChIP-seq, ChOR-seq, SCAR-seq, OK-seq, ATAC-seq, CUT&RUN, and CUT&Tag, and derivative protocols. CREPAS provides an end-to-end solution, from quality control to advanced analyses, including downsampling, peak calling, annotation, and visualization. By harnessing quantitative assays such as qChIP-seq and qChOR-seq, the normalization methods in CREPAS allow to compare the restoration kinetics of individual marks or proteins across replication timepoints. Moreover, the pipeline includes calculations such as fork directionality and partitioning using OK-seq and SCAR-seq data, linking replication dynamics to epigenetic inheritance. CREPAS is a valuable resource that enhances the efficiency and reproducibility of nascent chromatin sequencing data analyses, enabling the study of chromatin replication and propagation of epigenetic states. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=80 SRC="FIGDIR/small/732899v1_ufig1.gif" ALT="Figure 1"> View larger version (24K): org.highwire.dtl.DTLVardef@29ad26org.highwire.dtl.DTLVardef@26b9acorg.highwire.dtl.DTLVardef@67dcb8org.highwire.dtl.DTLVardef@cbc299_HPS_FORMAT_FIGEXP M_FIG C_FIG

Heterochromatin protein 1 (HP1 is an intrinsic component of heterochromatin domains where it is involved in a diverse set of functions including heterochromatin spreading and organization, chromatin compaction and transcriptional silencing. It has been suggested that HP1 functions through a phase separation mechanism, a process that has been observed in vitro in the presence of N-terminal phosphorylation, nucleic acids and nucleosome arrays. HP1 can also interact with numerous binding partners that contain a specific motif called an HP1 access code (HAC). HACs recognize and bind to an interface formed by the chromoshadow (CSD) domains in the HP1 homodimer, the functional form of the protein. It has been shown that some HP1 binding partners can enhance its phase separation ability while others disrupt the process. Here, we focus on the interactions between HP1 and three binding partners, namely the p150 subunit of the chromatin assembly factor 1 (CAF-1), the N-terminal domain of the lamin B receptor (LBR), and the mitotic protein Shugoshin 1 (Sgo1). Using phase separation assays, we show that CAF-1 prevents HP1 phase separation while LBR and Sgo1 enhance it. Binding assays, mutational studies, NMR spectroscopy and computational analysis allow us to dissect the contributions of the HAC motifs, the charge patterns of the binding partner sequences and the role of N-terminal phosphorylation on HP1 in condensate formation. Our results demonstrate that each binding partner uniquely balances these contributions to modulate the properties of HP1, while electrostatic interactions dominate the regulation of phosphorylated HP1. These results suggest that HP1 binding partners play an important role in the modulation of its properties and the regulation of its functions in distinct biological contexts.

Per- and polyfluoroalkyl substances (PFAS) are persistent toxicants with immunological, metabolic and epithelial effects, but their relevance to inflammatory skin disease remains unclear. We developed a computational toxicology framework to test whether perfluoroalkyl sulfonate programs, especially perfluorooctanesulfonic acid (PFOS) and perfluorohexanesulfonic acid (PFHxS), converge with psoriasis-associated keratinocyte inflammation. Exposure transcriptomes were derived from GSE236956, in which human embryonic stem cell-derived epithelial-lineage models were exposed to 10 M PFAS for 8-16 days. Six PFAS were prioritized using descriptors, Tanimoto similarity, toxicology evidence, adverse outcome pathway (AOP)-like key events, exposure differentially expressed gene burden and read-across support. PFAS signatures were integrated with psoriasis bulk transcriptomes, single-cell RNA sequencing, keratinocyte-state mapping, regulator and communication inference, spatial transcriptomics and token-aware Geneformer-compatible virtual perturbation. PFOS ranked highest in integrated prioritization, followed by PFHxS and perfluorooctanoic acid. PFHxS produced a smaller but directionally informative signature within a PFOS-dominant perfluoroalkyl sulfonate footprint. The shared PFOS and PFHxS program converged with psoriasis through inflammatory keratinocyte, epidermal-stress, cytoskeletal and lipid-related modules. Single-cell and spatial analyses localized the program to activated keratinocytes and inflammatory epidermal niches, with strong spatial co-localization with inflammatory keratinocyte and epidermal stress scores. Virtual perturbation prioritized S100A9, S100A8, KRT16, IL36G, CCL20, CXCL8, FABP5, KRT17, FOS, JUN and NFKBIZ as candidate effectors. These findings support an exposure-informed, experimentally testable hypothesis linking persistent perfluoroalkyl sulfonate programs to keratinocyte inflammatory niches in psoriasis.

In trypanosmatids genes, transcribed by RNA polymerase II do not have canonical promoters and are organized into directional gene clusters that mature into monocistronic transcripts by a co-transcriptional process known as trans-splicing. Even though gene expression is regulated mainly post-transcriptionally, it is currently understood that chromatin and epigenetics are also involved in this regulation. In eukaryotes, specific signals are normally required for the occurrence of an appropriate transcription initiation. Among them, trimethylation of histone H3 in lysine 4 is the most conserved signal normally detected at transcription start sites of actively transcribed genes. Unlike many model organisms, trypanosomes do not have defined promoters. Instead, transcription initiates in a bidirectional manner from dispersed regions coincident with divergent strand switch regions located between directional gene clusters (DGCs). In T. cruzi, H3K4me3 was observed at the origins of transcription coincident with divergent strand switch regions (dSSRs) in epimastigotes, but it has not been mapped throughout the whole genome at base-pair resolution or in other life stages so far. Here, we set up the CUT&RUN technique for T. cruzi epimastigotes and trypomastigotes. Consistent with a predominant post-transcriptional regulation along the life cycle, we did not find significant differences between life stages. We corroborated that H3K4me3 is enriched at dSSR adjacent to actively expressed DGCs. Moreover, we noticed that this histone mark exhibits different patterns that correlate with the genomic span of the transcription initiation regions and with transcriptional activity. Furthermore, we unveiled that the most actively transcribed DGCs are associated with shorter dSSRs and are located within the core compartment of the genome displaying a more accessible chromatin.

Interest in microbiota-host co-metabolism and the effects of its derived co-metabolites on biological processes is increasing rapidly. In addition to their demonstrated associations with mammalian metabolic health and cognition, microbiota-host co-metabolites (MHCMs) represent lifelong contributors to the endogenous exposome. We have previously shown the MHCM trimethylamine N-oxide (TMAO) to exert beneficial effects on murine blood-brain barrier integrity and cognition. Here we investigated whether these positive neural effects of TMAO extended to humans, analysing how TMAO exposure associates with neurodevelopmental outcomes in children and whether an in vitro human neuronal-astrocyte co-culture could contribute to further investigation of the underlying mechanism(s) and neuronal processes related to these associations. In a cohort study of childhood mental health (N=1,203), TMAO was associated with fewer internalising problems, while its precursor microbial metabolite trimethylamine was associated with more behavioural problems in both the cross-sectional and an independent longitudinal study from 1 to 15 years of age (N=630-820). Given prior associations between TMAO exposure and exposure to the environmental pollutants mercury and arsenic, we investigated how the effects of TMAO interacted with these known neurotoxicants. TMAO had a protective effect, modifying the relationship between arsenic exposure and poorer neurodevelopmental outcomes. Furthermore, TMAO activated synaptogenesis-related gene expression and was functionally protective against the negative effects of mercury in our in vitro model. Together, our findings emphasise the importance of interdisciplinary approaches to evaluate associations and potential pathways of MHCMs (endogenous) and environmental (exogenous) metabolites on neurodevelopment in exposome studies.

Histone proteins, which reside in the nuclei of eukaryotic cells, are involved in diverse cellular processes. The core histone H4 serves as a structural component of the nucleosome. Patients carrying mutations in H4Clustered histone (H4C) genes exhibit a broad spectrum of developmental abnormalities, including short stature, microcephaly, intellectual disability, growth retardation, and digital anomalies. However, the impact of H4 mutations on mammalian embryogenesis remains largely unclear. Here, we demonstrate that histone H4C genes play crucial roles in skeletal development and cortical neurogenesis. We found that mRNAs of the histone H4C gene family are specifically expressed in proliferating progenitor cells in the developing mouse neocortex and in human induced pluripotent stem cell-derived cortical organoids. CRISPR-mediated disruption of H4C3 in mice caused severe defects in skeletal formation and neocortical neurogenesis. Furthermore, overexpression of a mutant form of H4C3 resulted in altered expression of genes associated with cellular migration and motility. Together, these findings suggest that histone H4 plays a critical role in regulating the balance between proliferation and differentiation during mammalian embryonic development, thereby explaining the broad spectrum of patient phenotypes.

MYC oncoprotein binding at promoters and enhancers influences RNA polymerase II (RNAPII)-driven gene expression. Numerous genes also bind MYC near their transcriptional end sites (TESs). This often allows direct promoter-TES contact via looping and further regulates total and 'read-through' transcription that extends beyond standard termination sites. We aimed here to better clarify the rules governing TES associated MYC and/or RNAPII binding cross-talk in human and murine cells. Using ChIPseq and RNAseq datasets from the ENCODE portal and elsewhere, MYC and RNAPII binding profiles were found to differ around TESs and transcriptional start sites (TSSs). Variations in E box flanking sequences likely accounted for the somewhat lower affinities of MYC for TES-associated sites. Motifs for numerous other transcription factors were also observed to cluster non-randomly and in close proximity to MYC and RNAPII binding site peak summits. On average, genes with TES-proximal MYC or RNAPII sites were more highly expressed than those without, although co-binding tended to be suppressive. Both normal and neoplastic proliferative stimuli altered the MYC and RNAPII binding patterns of many genes, indicating that 'category switching' was common, subject to disparate external signals and often reversible. Functionally related gene sets with high levels of read-through transcription were uniformly marked by significant amounts of TES-associated MYC and/or RNAPII binding. These findings indicate that, both independently and together, MYC and RNAPII binding near TESs dynamically impact total and read-through transcription while also coordinating the expression of many common purpose gene sets.

SMAD4 is a central transcriptional effector of the TGF-{beta} signaling pathway and a frequently inactivated tumor suppressor gene in various cancers. Missense mutations in its MH2 domain are among the most prevalent somatic alterations in colorectal cancer (CRC). These mutations are associated with disease progression and poor prognosis, yet their precise mechanistic consequences have remained incompletely characterized. Here, we show that CRC-derived SMAD4 MH2 hotspot mutations (D351H, S357P, R361C, and R361H) selectively impair co-activator recruitment without disrupting chromatin occupancy. RNA-seq profiling demonstrated broad suppression of TGF-{beta} target gene expression across all mutants. Notably, the mutations confer distinct degrees of TGF-{beta} pathway unresponsiveness: R361H is completely refractory to TGF-{beta} stimulation, whereas R361C and S357P retain partial transcriptional responsiveness suggesting allele-specific differences in the severity of co-activator interface disruption. Genome-wide chromatin binding analysis by greenCUT&RUN confirmed that all mutants maintain wild-type-like genomic occupancy, as expected given that the MH1 DNA-binding domain is intact in each case. Proximity-dependent biotinylation mass spectrometry in COLO205 cells revealed that all four mutants exhibit markedly reduced interactions with the CREBBP/EP300 histone acetyltransferase complex and BRD4 relative to wild-type SMAD4 identifying disrupted co-activator engagement. Collectively, our findings establish that SMAD4 MH2 mutations impair TGF-{beta}-induced transcription by selectively reducing CREBBP/EP300 recruitment, which provides a molecular mechanism for the loss-of-function SMAD4 phenotype in CRC. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=154 SRC="FIGDIR/small/735541v1_ufig1.gif" ALT="Figure 1000"> View larger version (24K): org.highwire.dtl.DTLVardef@14f542eorg.highwire.dtl.DTLVardef@11fd220org.highwire.dtl.DTLVardef@1c3aa1org.highwire.dtl.DTLVardef@14d5a8e_HPS_FORMAT_FIGEXP M_FIG C_FIG

The sequence of the human genome provides a foundation for understanding cellular processes in health and disease. The organisation of this primary genetic information into cell-specific structure and function is critical to understanding the cell type-specific interpretation and execution of the genome. Epigenetic processes are essential for packaging and higher-level functional organisation of the genome, and changes therein are increasingly recognised as contributors to human disease. Building on primary data generated by multinational consortia, the International Human Epigenome Consortium (IHEC) has uniformly processed a collection of more than 2000 comprehensive human reference epigenomes, collectively referred to as EpiATLAS. This effort involved the development of standardised molecular and bioinformatics protocols, metadata models, and analytical tools to manage, integrate, display, and share vast amounts of epigenomic data. This includes the creation of a publicly available Epigenome Reference Registry, which provides a system for accessing protected human subject datasets and facilitates open searching of de-identified samples and experimental data. The integrated EpiATLAS ecosystem and its comprehensive human reference epigenome maps provide an unprecedented resource for the biosciences, expanding the annotated epigenomic landscape while uncovering previously unappreciated relationships among regulatory layers and revealing how epigenetic inputs underpin fundamental cellular functions and disease associations.

Blood-based DNA methylation may help identify biological changes related to dementia risk before clinical symptoms appear. We conducted an epigenome-wide association study of incident all-cause dementia in 5,999 cognitively healthy women from the Women's Health Initiative Memory Study, 777 of whom developed dementia over up to 25 years of follow-up. Baseline blood DNA methylation was tested for association with time to dementia. One CpG site, cg05917797, was significantly associated with dementia risk, with higher methylation linked to lower risk. This association was only minimally changed after accounting for APOE {varepsilon}4 carrier status, plasma p-tau217, and epigenetic aging measures. cg05917797 also replicated in meta-analysis of four independent prospective cohorts including 10,916 participants and 413 incident dementia cases. In post-mortem brain methylation datasets, higher methylation at cg05917797 was associated with lower Braak stage in temporal gyrus and cerebellum. In meta-analysis of all five prospective cohorts, including 16,915 participants and 1,190 dementia cases, cg05917797 remained the leading association, and three additional CpGs were identified for further study. These findings support cg05917797 as a reproducible blood-based epigenetic marker of long-term dementia risk.

IDH1 R132 mutations are among the most frequent hotspot mutations in cancer, but their mutational origins have remained unclear. Here, we provide evidence that IDH1 R132C, the predominant IDH1 mutation in cholangiocarcinoma, acute myeloid leukemia, and melanoma, likely arises through APOBEC3-mediated mutagenesis. IDH1 R132C is a TpC>TpT substitution on the lagging-strand DNA template within a hairpin-forming sequence context, consistent with APOBEC3 susceptibility. In vitro assays showed that APOBEC3A can deaminate the relevant cytosine, and APOBEC3A and APOBEC3B were relatively highly expressed in tumor types with recurrent IDH1 R132C mutations. IDH1 R132G, a TpC>TpG substitution at the same site, may similarly result from APOBEC3 activity. By contrast, IDH1 R132H, the predominant IDH1 mutation in lower grade glioma and glioblastoma, is a CpG>TpG substitution at a methylated cytosine on the leading-strand DNA template, a pattern more consistent with DNA polymerase epsilon replication error. Concordantly, tumor types enriched for IDH1 R132H showed relatively low POLE expression. Together, these in vitro and bioinformatic analyses provide insight into the distinct mutational mechanisms that likely underlie recurrent IDH1 hotspot mutations in cancer.

Epigenetic regulation is central to the developmental progression and pathogenicity of the unicellular eukaryotic parasite Plasmodium falciparum; yet, the contribution of DNA base modifications remains poorly understood. One such modification, 8-oxoguanine (8-oxoG), which was initially identified as an oxidative lesion and a marker of DNA damage, has since emerged as a transcriptional regulator in advanced eukaryotes. Given that P. falciparum encounters a highly oxidative environment in human blood, we investigated the potential gene regulatory role of 8-oxoG during its intra-erythrocytic developmental cycle (IDC). Using immunodetection assays, we first confirmed the presence of 8-oxoG in P. falciparum genomic DNA and observed a gradual increase in 8-oxoG abundance from ring to schizont stages. We then optimized oxidative DNA immunoprecipitation sequencing (OxiDIP-seq) for the highly AT-rich parasite genome and generated genome-wide 8-oxoG profiles across four IDC timepoints, which revealed reproducible enrichment of 8-oxoG at discrete genomic loci, with more than 50% of the peaks stable across developmental stages. Notably, 8-oxoG accumulated at putative G-quadruplex-forming sequences in the parasite genome and preferentially localized within exonic regions of protein-coding genes, exhibiting a marked enrichment near STOP codons and within 3' untranslated regions. This in turn correlated with significantly higher steady-state transcript levels of 8-oxoG-marked genes, with stage-specific changes in 8-oxoG enrichment closely matching transcriptional activity. Furthermore, 8-oxoG-marked loci were preferentially associated with active and poised histone post-translational modifications, while showing no evidence of altered nucleosome occupancy. Collectively, these findings demonstrate that 8-oxoG is a widespread and non-random DNA modification in P. falciparum and suggest that it may function as an epigenetic mark associated with transcriptionally permissive chromatin and gene activation during parasite blood-stage development.

Tenuazonic acid (TeA) is an emerging Alternaria mycotoxin frequently detected in food and feed commodities, raising concerns about its toxicological relevance. Chronic oral exposure to TeA has been reported to induce dysplastic alterations in the esophageal mucosa of mice, while human biomonitoring data indicate an association between TeA exposure and esophageal cancer, although a causal relationship has not yet been established. At a mechanistic level, the effects of TeA in esophageal cells remain poorly characterized. Therefore, this study investigated the impact of TeA on cytotoxicity, oxidative stress, DNA damage, mitochondrial homeostasis, cell-cycle distribution and transcriptomic stress responses in human esophageal KYSE-510 cells. TeA induced a concentration-dependent reduction in metabolic activity and total protein content after 24 h exposure to 0.1-100 M. Significant cytotoxicity was measured starting from 20 M. At sub-cytotoxic concentrations, TeA triggered rapid ROS formation within 5-30 min exposure and induced formamidopyrimidine-DNA glycosylase (FPG) sensitive DNA damage after 1 h exposure (5-7.5 M), indicating oxidative DNA lesions. In addition, TeA altered mitochondrial morphology after 4 h exposure at 7.5 M, manifested by shrinkage of the mitochondrial network area and perinuclear redistribution, while mitochondrial respiration showed only a non-significant tendency towards reduced respiratory capacity. RNA sequencing after 6 h exposure to 10 M TeA revealed oxidative stress-associated transcriptional changes, impaired antioxidant and stress-adaptive responses, and p53-associated stress signaling. Furthermore, TeA induced significant G2/M phase accumulation after 24 h exposure to 1-10 M.

TRIM28 is an epigenetic co-repressor protein that silences transposable elements (TEs). Although loss-of-function studies in mice led to neurodevelopmental defects, a functional link between TRIM28 and human neurodevelopment has yet to be established. In this study, we describe two patients with neurodevelopmental delay who carry de novo TRIM28 missense variants. Using CRISPR-edited induced pluripotent stem cell lines and differentiated neural organoids, we demonstrate that these variants result in the loss of the histone mark H3K9me3 over TEs. This releases the regulatory potential of TEs resulting in altered expression of nearby genes. These findings could be replicated using CRISPRi-based TRIM28 silencing, which suggests that the two variants result in a loss of function. Our results highlight the critical role of TRIM28 in regulating TEs during human brain development, establishing a link between TRIM28 variants and neurodevelopmental delay.

Mammalian sperm chromatin carries epigenetic information with the potential to influence offspring phenotype, making its faithful characterization essential. It has been suggested that cauda sperm preparations are contaminated by somatic chromatin, that this contamination dominates genome-wide profiles, and that valid results require pretreatment with somatic cell lysis buffer, DNase I, and dithiothreitol. Here we show that properly purified cauda sperm contain no detectable somatic cells or cell-free DNA and that this pretreatment disrupts sperm chromatin organization. SCLB permeabilizes the sperm nucleus, allowing DNase I to fragment the sperm genome in situ, while DTT treatment causes chromatin to leak out of the nucleus. Using ATAC-see, we further demonstrate that Tn5 transposase can access intact protamine-condensed sperm chromatin without DTT, refuting the premise that profiles from untreated sperm reflect contamination. Pretreatment therefore damages the chromatin it claims to purify, and published profiles of untreated cauda sperm are valid and require no systematic re-examination.

Background Complex traits are influenced by numerous variants, most of which have regulatory effects on gene expression that can be mediated by DNA methylation. Molecular QTL mapping is an approach that aims to dissect these effects. However, obtaining molecular phenotypes on a large scale is challenging, particularly in livestock species. In cattle, an epigenotyping array called EpiChip has recently been developed in the European RUMIGEN project. The EpiChip, which contains 43,317 CpG sites distributed all over the bovine genome, enables large-scale measurement of DNA methylation. This study aims to characterize the genetic determinism of blood DNA methylation in cows by estimating heritability and mapping cis- and trans-methylation QTLs (meQTLs). Results Whole blood samples from 4,457 genotyped Holstein cows were epigenotyped. Across all CpG sites, the heritability estimates averaged 24.6%. The local meQTL mapping at sequence-level for variable CpG sites (SD > 2.5%; n = 28,806) detected cis-meQTLs for 80.1% of the CpG sites, with sentinel SNPs located close to their associated CpGs. A two-step analysis was also conducted to identify long-range associations, with a particular focus on trans-meQTL hotspots. First, we identified CpG-SNP trans-associations using medium-density genotypes (50k SNPs) that revealed 31,846 SNPs with significant effects on 1 to 530 trans-CpG sites. Then, regions associated with at least 34 independent trans-CpGs were retained defining 31 hotpots. For each hotspot, a local sequence-level GWAS was conducted using the first principal component derived from the associated trans-CpGs. Out of the 31 detected hotspots, three were located close to transcription factor genes (RUNX1, NFIC and FOXA3) for which the associated trans-CpGs were enriched for the corresponding binding motif. Two other hotspots were located within KDM5A and KDM5B, and their corresponding trans-CpGs were strongly overrepresented in H3K4me3 narrow peaks in blood as well as in other tissues. Conclusions By identifying functional candidate genes associated with blood DNA methylation in cattle, these findings provide new insights into the regulatory architecture of DNA methylation in mammals, highlighting the value of large-scale molecular data from livestock populations.